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We theoretically study the joint influence of uniform Dzyaloshinskii-Moriya (DM) interactions, 
symmetric exchange anisotropy (with its axis parallel to the DM vector) and arbitrarily oriented 
magnetic fields on one-dimensional spin 1/2 antiferromagnets. We show that the zero-temperature 
phase diagram contains three competing phases: (i) an antiferromagnet with Neel vector in the plane 
spanned by the DM vector and the magnetic field, (ii) a dimerized antiferromagnet with Neel vector 
perpendicular to both the DM vector and the magnetic field, and (iii) a gapless Luttinger liquid. 
Phase (i) is destroyed by a small magnetic field component along the DM vector and is furthermore 
unstable beyond a critical value of easy-plane anisotropy, which we estimate using Abelian and non- 
Abelian bosonization along with perturbative renormalization group. We propose a mathematical 
equivalent of the spin model in a one-dimensional Josephson junction (JJ) array located in proximity 
to a bulk superconductor. We discuss the analogues of the magnetic phases in the superconducting 
context and comment on their experimental viability. 
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I. INTRODUCTION 

Quantum spin chains continue to be an active focus of 
research largely because they serve as interesting model 
systems to explore strongly correlated quantum order 
in low dimensional antiferromagnets^ superconductors^ 
and ultracold atoms.— A significant fraction of current 
research is devoted to the understanding of frustrated 
quantum magnets, which display a host of exotic ground 
states^ One of the agents responsible for magnetic frus- 
tration is the Dzyaloshinskii-Moriya^ interaction, • 
(Si x Sj), which originates from spin-orbit coupling and 
broken inversion symmetry is the spin operator at site 
i and Djj is the DM vector. Since this interaction may 
induce spiral spin arrangements in the ground stated it 
is intertwined with ferroelectricity in multiferroic spin 
chains^ Besides, the DM interaction plays an important 
role in explaining the electron spin resonance experiments 
in some one-dimensional antiferromagnets^ Moreover, 
the DM interaction modifies the dynamic properties^ 
and quantum entanglement^ of spin chains. 

The present work is motivated by an elegant recent 
studyi^ that has predicted intriguing field-induced anti- 
ferromagnetic order in Heisenberg spin 1/2 chains with 
uniform DM interaction (Dij = Dz for any A cru- 
cial aspect of Ref. [l2[ is that the external magnetic field 
is taken to be (nearly) transverse to the DM vector, so 
that the system has fully broken spin rotational symme- 
try. The main objectives of our work are to generalize 
the analysis of Ref. [l2j for the case of nonzero symmet- 
ric exchange anisotropy as well as to find new physical 
contexts where it might be experimentally testable. 

We begin in Section II by introducing the pertinent 
spin model. In Section III we identify the classical 
ground states using the large-spin approximation. Al- 



though some features predicted by the classical analysis 
are erroneous, most of the classical ground states have 
a correspondent in the more rigorous quantum analysis 
performed in Section IV. In particular, we derive a new 
result regarding how a magnetic field component along 
the DM vector modifies the classical soliton-lattice phase. 

Section IV is divided in two subsections which com- 
plement each other to an extent. The first subsection 
approaches the problem from a non-Abelian bosoniza- 
tion perspective and is constructed around the elegant 
chiral rotation introduced in Refs. [l2llT3j . This sub- 
section constitutes the core of the present work. The 
second subsection revisits the problem from an Abelian 
bosonization viewpoint, and may be skipped on a first 
reading without loss of continuity. The phase diagram is 
richest in the case of weak easy-plane anisotropy. In this 
regime there are three competing ground states: (i) an 
antiferromagnet with Neel vector in the plane spanned 
by the DM vector and the magnetic field, (ii) a dimer- 
ized antiferromagnet with Neel vector perpendicular to 
both the DM vector and the magnetic field, and (iii) a 
gapless Luttinger liquid. Phase (i) was first identified 
in Ref. [l2| at the isotropic exchange point. We demon- 
strate that this phase is fragile under weak-to-moderate 
easy-plane symmetric exchange anisotropy, and estimate 
the critical value of the anisotropy beyond which it disap- 
pears. Furthermore, we show that a very small magnetic 
field component along the direction of the DM vector 
suffices to destabilize phase (i) in the neighborhood of 
the isotropic exchange point. Thus the experimental de- 
tectability of phase (i) appears unlikely except in chains 
with easy-axis anisotropy. Phases (ii) and (iii) emerge 
as a consequence of symmetric exchange anisotropy and 
constitute the main findings of this work. The Luttinger 
liquid ground state prevails when the DM interaction is 
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large compared to the magnetic field perpendicular to 
the DM vector, and the antiferromagnetic phase is sta- 
bilized under the opposite condition. Even though the 
outcomes of the Abelian and non-Abelian methods agree 
roughly, the former method misses a few key features 
such as the coexistence of antifcrromagnetism and dimer- 
ization in phase (ii). This is more a merit of the chiral 
rotation carried out in conjunction with the non-Abelian 
treatment than an intrinsic flaw of the Abelian bosoniza- 
tion. 

In Section V we recast the magnetic model onto 
a mathematically equivalent problem that consists of 
a one-dimensional Josephson junction array located in 
close proximity to a bulk superconductor and placed un- 
der a magnetic field. The analysis of Section IV can be di- 
rectly transferred to determine the phase diagram of this 
system at small or large magnetic fields, depending on 
whether the array is made of 7r-junctions or conventional 
junctions, respectively. Phases (i) and (ii) discussed in 
Section IV correspond to a a charge density wave and a 
vortex lattice, respectively. In the vortex lattice phase, 
circulating supercurrents flow between the array and the 
bulk superconductor. The magnitude of these circulating 
currents oscillates from one "plaquette" to another; this 
is how dimerization manifests itself in the superconduct- 
ing context. The magnetic Luttinger liquid maps into a 
state in which the Josephson coupling between the super- 
conducting islands and the bulk superconductor becomes 
irrelevant. The transitions between these phases may be 
controlled with external magnetic fields and by engineer- 
ing material parameters. We outline the desiderata for 
an experimental implementation. 

Section VI contains a brief summary of this work and 
the Appendices include a few technical details concerning 
the classical and quantum phase diagrams. 

II. MODEL 

Consider an S = 1/2 one-dimensional antiferromag- 
netic chain in presence of a uniform Dzyaloshinskii- 
Moriya (DM) interaction and an external magnetic field. 
Its Hamiltonian is 

i 

+Dz ■ (Sj x S j+1 ) - Y,( h * SX + h ^S z ), (2.1) 

3 3 

where J is the exchange coupling, A is the symmetric 
exchange anisotropy parameter, D is the strength of the 
DM interaction, z is the direction of the DM vector (cho- 
sen to be parallel to the symmetric exchange anisotropy 
axis), and h x (h z ) is the component of the magnetic field 
perpendicular (parallel) to the DM vector. This model 
describes one-dimensional magnetic systems with broken 
inversion symmetry, as well as interacting quantum wires 
with spin-orbit interactions.— As we demonstrate in Sec- 
tion V, Eq. (|2.1[) is also germane for a one-dimensional 



array of Josephson junctions that are proximity coupled 
to a bulk superconductor (Fig. [T5|). 

Eq. (|2.1[) may be rewritten in a physically more sugges- 
tive manner by rotating the spins as Sj = exp(— iaj)S^~ , 
where a = tan -1 (D/J) S [0,7r/2]. This rotation gauges 
away the DM interaction and produces a XXZ antiferro- 
magnet with an altered exchange anisotropy A c fj and a 
magnetic field that rotates with a pitch angle a in the 
plane perpendicular to the DM vector: 

3 

- h * E (• s "; , " w + ~sj^ ia3 ) - h z y, Sj, ( 2 - 2 ) 

3 3 

where J = J/ cos a and A c jf = A cos a < A. When 
D > J the spiral magnetic fields rotates rapidly and thus 
Eq. (|2.2[) can be mapped onto a XXZ model with a mag- 
netic field along z. When D ~ the spiral magnetic 
field rotates very slowly and Eq. (|A8|) transforms onto 
a XXZ model with a spatially uniform magnetic field in 
the xz plane. In this work we shall be concerned with 
< D,h x ,h z « J. 

III. CLASSICAL ANALYSIS 

Eq. (|2.ip is a complicated model with fully broken spin 
rotational symmetry. For pedagogical purposes it is use- 
ful to begin with simple classical considerations which 
shed light on the possible ground states of the fully quan- 
tum mechanical problem. When S is large, it is adequate 
to substitute Sj = S(sin9j cos <f>j, sin 9 3 ; sin 4>j, cos 9j) in 
Eq. (|2.1|) and seek solutions that satisfy dH/dfa = 
dTi/dOi = for all i. Wc limit ourselves to A e s > 0. 
For h z = we find five distinct phases: 

(i) Uniform ferromagnet (FM), i.e. fa = (aligned 
with the field) and 0j = 9. There are two solutions: 
9 = tt/2 or 9 = sm~ 1 (h x /2JS(l - A)). Its energy per 
site is e FM = JS 2 [A + (1 - A) sin 2 9} - h x S sin 9. 

(ii) Uniform antifcrromagnet with Neel vector along 
x ("N x "), i.e. 9j = tt/2 and <fij = irj. Its energy per 
particle is 6n x = —JS 2 . 

(iii) Uniform antiferromagnet with Neel vector along 
y ("JV»"), i.e. 0j = tt/2 and fa = {-1)^^/2 - <j) Q ). 
4>q = sin -1 (h x I '4J S) is the canting angle towards the 
direction of the magnetic field (x). Its energy per site is 
£Ny = —JS 2 — h 2 /8J. 6N y < £nx whenever h x ^ 0. 

(iv) Uniform antiferromagnet with Neel vector along 
£ ("iV z "), i.e. Oj = (— 1) J '6> + irj and fa = 0. 9 = 
sin~ 1 (h x /2JS(l + A)) is the canting angle towards the 
direction of the magnetic field. The energy per site for 
this phase is e Nz = -JS 2 A - /i 2 /4J(A + 1). 

(v) Spiral XY antiferromagnet (LL). For h x = 0, 
this phase is characterized by 9j = tt/2 and fa — 
a .i + + X- X is the- global angle of the spiral; its 
arbitrariness renders the LL phase gapless, in contrast 
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FIG. 1: Classical magnetic configuration in the "iV 2 "' (a) and 
LL (b) phases when h x 7^ 0, D 7^ and h z = 0. All spins 
lie in the a;y plane. The solid arrows replicate the staggered 
part of the magnetization whereas the dotted arrows represent 
the uniform canting towards the direction of the field. The 
uniform component of the spins is spatially inhomogeneous 
in the LL phase and it vanishes at the center of the solitons. 
In (b) we limit ourselves to one soliton; for < h x < nDS 
these solitons form a periodic array^ with a periodicity that 
is incommensurate with the underlying lattice. 



to the ones introduced above. Its energy per site is 
6ll = —JS 2 /cosa. For h x 7^ 0, the spiral distorts into 
an incommensurate soliton latticed with 9j = tt/2. A 
single soliton is described by <pj = nj + a(j) + (— 
where a(j) = 2tan _1 [exp(j/i a; /2J5')] — tt/2 and b(j) — 
— h x /2JS tanh(j7i x /2 JS) vary slowly on the scale of a 
lattice spacing. Fig. Q] illustrates the spin arrangement in 
the soliton. 

When h x = 0, it is easy to verify analytically that the 
classical ground state is LL (if A c ff = A cos a < 1) or N z 
(if A eff > 1). When h x e (Q,irDS) and A < 1 the classi- 
cal ground state is LL (incommensurate soliton lattice). 
h x = ttDS is the critical field for the commensurate- 
incommensurate transition. This critical field is inde- 
pendent of the strength of easy-plane anisotropy because 
0j = n/2. For h x € (nDS, AS J) and A < 1 the ground 
state is "iV 3 "' . When h x > AS J the ground state is FM. 
For the purposes of this paper h x will never be large 
enough to stabilize the FM phase (but see Section V 
for an exception). For A > 1, the ground state is "N z " 
regardless h x . 

Thus far we have neglected h z . After turning on 
h z , " N v " develops a uniform canting towards z with 
cos6> ~ h z /2JS(l + A). On the other hand the soli- 
ton develops a canting that has both uniform and stag- 
gered components. As we explain in Appendix A, this 
leads to a redefinition of the soliton parameters that re- 
sults in an increased critical field for the commensurate- 
incommensurate transition: 



h x . 



tJS 



8J^2 (1 _ A 2 ff) 



(3.1) 



Eq. d3U applies when |A-1| >> h 2 J J 2 S 2 , (da/dx-a) 2 . 
It also indicates that the influence of h z on the critical 
field gets weaker when A decreases. Regarding l W z ," 
its Neel vector tilts towards x by an angle 0. For A ~ 
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FIG. 2: Classical magnetic configuration when h x 7^ 0, D 7^ 
and h z 7^ 0. (a) "jV y ": the dotted circles indicate uniform 
canting towards z. (b) LL: the spins conform into a soli- 
ton lattice in the xy plane much like in Fig. [T] but they are 
now canted towards z as well. As detailed in Appendix A, 
the latter canting is spatially inhomogeneous. On one hand, 
it has a non-staggered component (not shown in this figure) 
that varies slowly along the soliton. In addition, the canting 
angle has a staggered component (represented in the figure 
via dotted and crossed circles) that is nonzero only at and 
near the core of the soliton. (c) " N z ": due to h z 7^ the Neel 
vector is tilted from z towards x. Consequently the canting 
component is not uniform. The crossed circles denote a mag- 
netization component that points into the page. 



1, sin 6 ~ —h z /y?h x + hi and the Neel vector is nearly 
perpendicular to h = h x x + h z z. Fig. [2] illustrates the 
influence of h z in the classical ground states. 

Turning on h z modifies the phase diagram qualita- 
tively. For A < 1, LL and U N V " remain as the ground 
states, although the critical field for the spin-flop transi- 
tion is now A-depcndcnt. For A > 1 the ground state is 
l W z " with a tilted Neel vector. 

Some of these classical considerations apply in the 
case of the 5 = 1/2 spin chains, while others do not. 
For instance, the occurrence of "N v " and gapless (spi- 
ral) ground states will be corroborated by the upcoming- 
quantum analysis. Similarly, our study will confirm that 
the critical field for the commensurate-incommensurate 
transition is independent of A at h z = but not at 
h z =/= 0. In contrast, the quantum analysis will show 
that "jV z " can be the ground state even at A < 1, thus 
refuting the classical prediction. 



IV. QUANTUM ANALYSIS 

In this section we analyze the exact low energy behav- 
ior of Eq. (|2.1j) by combining rcnormalization group argu- 
ments with Abelian and non-Abelian bosonization tech- 
niques. The former technique is best suited for large easy 
plane anisotropics (A e g ~ 0) whereas the latter method 
is most reliable at A e g = 1. Bearing in mind that each 
approach has its shortcomings, we shall compare them 
with each other when possible. 
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A. Non-Abelian bosonization 

Many S = 1/2 antifcrromagnetic chains showcase a 
weak exchange anisotropy because they involve Cu 2+ 
compounds. Furthermore, spin-orbit coupled antifcr- 
romagnets without inversion symmetry typically ex- 
hibit A e ff ~ 1 because the anisotropy induced by the 
DM interaction nearly cancels the preexisting exchange 
anisotropy*^ 

The ground state properties of these nearly isotropic 
systems may be conveniently accessed using the non- 
Abelian bosonization ^ 16 ' 17 In this framework, the spin 
operators are approximated as 



J L (x)+3 R (x) + (-l) x ' a N(x) 



(4.1) 



where a is the lattice constant. J/, and 3 R are the uni- 
form components of the left- (right-) spin-currents, re- 
spectively. These currents are conserved at the SU(2) 
symmetric point (A c ff = 1). N is the staggered compo- 
nent of the local spin density. 

In the continuum limit Eq. (|2.1[) is bosonized in terms 
of the SU{2)\ Wess-Zumino-Novikov-Witten model. The 
low energy effective Hamiltonian can be written in the 
Sugawara form, which is quadratic in the SU(2) currents: 



H = Ho+ Hbs + V, 



(4.2) 



where 
"Ho 



R ' J R 



■Jr. - J; 



J dx[r R r L + jpi + (i + \ xc )r R J z L } 



V = -h x dx{J x R + Jt) - K \ dx(J R + J 



+D J dx{J z R - Jf ) - 3bs A DM J dxjyUA.Z) 

where v ~ Jair/2 is the velocity of the low energy excita- 
tions near the isotropic point (we neglect the anisotropy- 
induced rcnormalization of the velocity). "Ho is the 
non-interacting part, the backscattering part "Hbs is the 
leading marginally irrelevant interaction, and V collects 
the second line of Eq. (|2.1[) . g hs is the (temperature- 
dependent) coupling constant for the effective interac- 
tion, A xc = c(l— A) is the symmetric exchange anisotropy 
parameter and c is a positive constant. Also, D = 
D(l + 2j 2 )/ir and Adm = dD 2 /J 2 , where d > is an- 
other constant. 7 ~ O(l) is the mean-field expectation 
value of the charge operator. 



A = A xc + A 



DM 



(4.4) 



is the effective anisotropy parameter, such that A < 
corresponds to easy-axis anisotropy (easy axis: z) and 



A > indicates an easy-plane anisotropy (easy plane: 
xy). 

The marginal coupling g^s at energy scale T was de- 
termined accurately from the Bethe ansatz^ 



^y + 5log[ 5 (T)] = lo, 



T 



(4.5) 



where g(T) = <7bs/27ru, 7 = 0.577... and T is the temper- 
ature. 

While deriving the expression for V in Eq. (j4.3[) we have 
exploited the operator product expansion^ 7 - (OPE) and 
have kept a higher order term (oc D 2 /J 2 ) when expand- 
ing in gradients of d x <3. This term has been neglected in 
previous studies^ yet as we shall show below it is essen- 
tial in order to reproduce the correct phase diagram in 
some simple limits. Hereafter we absorb this term into 

n hs . 

The constants c,d and 7 are related to each other 
and may be determined exactly in certain limits. For 
instance, we know that the total anisotropy parameter 
A = X X c + Adm must vanish when A e g = 1, i.e. when 
A = (1 + D 2 / J 2 ) 1 / 2 . Expanding the square root to lead- 
ing order in D/J, we arrive at d = c/2. Furthermore, 
when h x = and A e ff < 1 it is well-known (e.g. from 
the Bethe ansatz method) that the ground state is a Lut- 
tingcr liquid with a specific value of the Luttinger param- 
eter. This determines c in absence of fields (more on this 
below). Finally, the value of 7 may be extracted via the 
OPE of the DM interaction, which relates 7 with d . For 
the purposes of this work, the precise value of 7 will not 
be important and we shall not be concerned about the 
distinction between D and D in our numerical calcula- 
tions. The renormalization group analysis of Eq. (|4.2[) 
can be simplified considerably by applying a chiral rota- 
tion around the y axis^, which acts differently on right 
and left currents: 



= K(6 R{L) )M R{L) 



cos 9 



R(L) 



sm9 R{L 



1 







= tan 



sin6» fl(L) cos9 R(L) 



D 



9l = 



tan 



1 ( 



d L = D + ft 2 (4.6) 



The key motivation for the chiral rotation is that it re- 
casts V into 



V = dx 



d 2 R + h 2 x M* R 



(4.7) 



which is an operator proportional to M 



r/l- 



5 



Consequently it can be eliminated^ by the following position-dependent phase shifts: 



M+ = Ml + iM\ 



M R exp [—i{t^ - t e )a 
AI+ = Ml + iM v L -> M+ exp + tg)x] 



M R -+ M R 



td 



to 



Ml 



Ml 



4-7T 
td, + tg 



4tt 



(4.8) 



where 



U = U*L + h l + \f%l + h i) /2« 



fe = 



+ hi - Jdl + hi ) /2i 



(4.9) 



The fact that the DM interaction and the magnetic field may be treated on the same footing and absorbed together 
through phase shifts is a qualitative advantage of the chiral rotation.— While Ho remains invariant throughout, the 
successive trasformations modify the backscattering term as follows: 



Tibs = 2nv / dx 
where 



y A {M R Ml-M%Ml)+y A {Ml i Ml + M R Ml)+ £ y a M%Ml 



-^■Ha + Ha + H! 



'Ha, 



n A 
h b 

He 



ttv y A J dx(M R M+e lit f +te)x - M^Af£e~^**~* 9 ^ + h.c.) 
ttv y A J dx(M R M+e i(t * +te)x + M+Mler l{t *-~ te)x + h.c.) 
ttv y B J dx(M+M£e- l2t t> x + h.c.) 
ttv y c [ dx(M+M+e l2teX + h.c.) 



-2ttv y a 



dxM R Ml 



(4.10) 



In Eq. (|4.10j) we have neglected small terms that originate from the shifts in M£ and M R . The initial values for the 
coupling constants in Eq. (|4.10p are given by 



lfe(0) 

y y (o) 
vM 

Va(P) 

va(o) 



9bs 
2ttv 

flbs 

2ttv 

ffbs 

2nv 



1 



A 



1 + - ) cos ( 



A 

--COS* 



-cos 6+ 
2 



9bs 
2ttv 



sin 8 



9bs A 

sin f 

2irv 2 



(4.11) 



where 6' ± = 9 R ± 9l and 



yc 

Vb 

Ha 



[Vx - Vy) 



-y z - 
r 



(4.12) 



The anisotropy parameter A modifies the initial couplings and combines with h z to introduce an extra coupling con- 
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stant yA in the RG equations (note that J/a(0) = when 
h, = 0). 

Because of the oscillatory phase factors introduced by 
Eq. (|4.8p . the RG analysis of H must be carried out in 
multiple stages. In the hrst stage we integrate out mo- 
menta that are large compared to max(t^,,tg) and thus 
all phase factors may be ignored in "Hbs- The flow equa- 
tions can then be derived in the standard manner— using 
OPE: 



dy x 
dl 

dy y 
dl 

dy z 
dl 

dyA 
dl 

dyA 
dl 



VzVy 

y z y x - (va + VA)(yA - va) 

VxVy 

VyVA 
-VyVA- 



(4.13) 



For convenience we begin integrating Eq. (|4. 13[) at an 
initial energy scale To = 0.077J, where the effective cou- 
pling has the value <7b s (?o) ~ 0.23 x (2itv) as dictated 
by Eq. (|4.5j) . This is a low enough energy scale, and 
a small enough value of <7bs(T), that the above (lowest 
order) RG equations apply, at least approximately. We 
then integrate Eq. (|4.13|) towards lower energy scales in 
order to determine the zero-temperature phase diagram. 
I = log(T/ao), where L is the length of the chain and 
ao = v/Tq = 20.4a is the ultraviolet RG cutoff length- 
scale. 

When A = Eqs. (j4~TTj) and (|4~T3)) reduce to those 
shown in Ref. (T^]. The fact that A > for an isotropic 
(A xc = 0) Heisenberg antiferromagnet appears to have 
been overlooked by Ref. [IH, which takes A xc = and 
A c ff = 1 simultaneously even in presence of DM interac- 
tions. 

Eq. (|4.13[) is no longer valid when I > 
min(log(l/aot0),log(l/aofe)). For definitcness we 
assume tg << which holds e.g. when h z << D,h x . 
Then, at I > 1$ = log(l/aoi0), exp(it<f,x) oscillates 
rapidly and thus the factors multiplying it in %bs aver- 
age to zero. Therefore yA, jjA and ys stop renormalizing 
at I = 1$. In contrast, the factor that multiplies yc 
is approximately uniform because tg << t$. The flow 
equations for the second RG stage are obtained by 
setti ng y A = Va = and y B = (i.e. y x = —y y ) in 
Eq. gT3J|. This yields 



dye 

dl 
dy a 

dl 



y<jyc 



y 2 c, 



(4.14) 



which are the famous Kosterlitz-Thouless (KT) equa- 
tions with known analytic solution. The "initial" condi- 
tions for the second RG stage are y a (l<j>) = —yz{l<j>) and 
yc{l<t>) = {Vx{l<l>) ~ y y {l<j>))/2. We integrate Eq. (|4.14j) 



from I = Ic/) up until I = lg = log(l/aofe). At I = lg 
the flow of yc stops because its coefficient in Tibs con- 
tains a exp(itgx) factor. By setting yc = in Eq. I|4.14|) . 
it follows that y a stops flowing as well. Thus at I = 1$ 
the system is in its ground state, the nature of which is 
determined by the final values of the coupling constants. 

In order to elicit the physical meaning of the different 
ground states we resort to the relation^ between the 
non-Abelian operators and the bosonic fields (0, <&) with 
which the Abelian bosonization is constructed: 



Ml 

Af ± 

N z 



2na ' L 2na 



1 



2V2tt 



(3 X $ + 3 X 6) ; M Z R 



1 



2V27T 



7 +i v / 2i : e 
7ra 

7 



TV (I 



sin(V27r<I>) 



(d x $ - d x e) 



(4.15) 



where (Af x ,J\f y ,J\f z ) is the staggered magnetization in 
the rotated frame. <& and 9^8 are canonically conjugate 
fields, i.e. [$>(x),d x Q(x')] — iS{x — x'), We reserve the 
and $ notation (without tilde) for the next subsection, 
where we shall employ Abelian bosonization in a differ- 
ent coordinate system. Eq. (|4. 1 5[) yields Mi Mf +h.c. oc 



cos(v / 87rG>) and M^M£ + h.c. cx cos(v / 87r ( i ) ). Hence if 
Uc(Je) ~^ ioo the minimum energy state corresponds to 
cos(-\/87rO) = =pl, which implies that 9 gets ordered and 
a gap is opened in the spin excitation spectrum. Further- 
more, Eq. (|4.15|) shows that when cos(v / 87T0) = 1(— 1) 
there is long-range antifcrromagnctic order with the Nccl 
vector pointing along x (y) in the rotated frame. 

When \yB{le)\ oo the minimum energy state cor- 
responds to cos(v / 87r'i > ) = T-l- Once again resorting to 
Eq. (|4.15j) . it follows that when cos(v / 87r$) = lor - 1 
the system settles into a dimerized state or an antifcr- 
romagnetic state with Neel vector along z in the rotated 
frame, respectively. Both phases are gapped because $ 
is ordered. 

In the LL phase neither O nor $ order. This gap- 
less Luttinger Liquid has dominant dimer or spin den- 
sity wave correlations depending on whether yc{le) dom- 
inates over yB(le) or vice versa. At first glance a gap- 
less ground state appears unlikely in a model with com- 
pletely broken spin rotational symmetry; however, it is 
not unprecedented. Similar states with gaplcss spin ex- 
citation spectra occur in magnetized, spin-orbit coupled, 
one-dimensional conductors that are placed under a mag- 
netic field.— As we shall see below, the LL ground state 
results when D dominates over h x . 

In order to determine what the aforementioned ordered 
states mean in terms of the original spin variables J and 
N, we once again perform the rotation introduced in Eq. 
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(14. 6|) and arrive at 

TV* 

AT 2 



' t) m * ~ sin (t 1 

[ Af x + COS f^— ) 



[ — J Af v - sin ^ — ) 



sin ( — J M y + cos f — ] £ . 



(4.16) 



J\f % and £ denote the Neel vectors and the dimerization in 
the rotated frame, respectively. N l and e are their coun- 
terparts in the original coordinate system. The deriva- 
tion of Eq. (|4.16[) becomes straightforward after recog- 
nizing tha1>^ 

M R oc tY{crg~ 1 d z g} ; M L oc tx[(rg~ l d s g] 
N = tr[ergr] 

£ = tr[g] (4.17) 

where a is a vector of Pauli matrices, d z ^) = 9 t /v+(—)d x 
and 



g oc 




FIG. 3: Phase diagram corresponding to Eq. (|2.ip : derived 
using non-Abelian bosonization. We fix D = 5 x 10 -4 J and 
h z = 0; we vary h x and A e ff (note that A e ff and A are nearly 
identical in this figure). For A c ff > 1 (easy axis anisotropy), 
the ground state is unequivocally "N z " . For A e ff < 1 (easy 
plane anisotropy), two new phases emerge: the gapless LL 
at D » h x and the gapped u N y " at D « h x . The phase 
line separating "N v " and "A z " has a horizontal asymptote 
at Aefr = 1 as h x /D — > oo. Below a critical value of the 
anisotropy (A c ff < A c ~ 0.973 in this figure) "A z " can no 
longer be the ground state. 



(4.18) 



is the SU(2) matrix field that enters the WZNW action 
and transforms as 



g -+ e^Wge 



i/2 



(4.19) 



under a chiral rotation (this g is of course not to be con- 
fused with its homonym of Eq. (|4.5|) ) . While we find 
(£} = for all values of D, h x and h z , there exist regions 
of parameter space for which (M y ) ^ 0. This then trans- 
lates into a coexistence of antiferromagnetism (with Neel 
vector along y) and dimerization in the original frame. 
Table 1 enumerates and characterizes this and other pos- 
sible ground states of Eq. (|2.1[) in the original coordinate 
system. 

In Appendix B we digress on some simple limits in 
which the phase diagram is known with certainty. Be- 
sides providing a reality check, this enables us to de- 
termine the value of the constant c introduced above 
Eq. (14.41) . This appendix may be skipped on a first read- 
ing. 

We now embark on the numerical study of the general 
phase diagram for Eq. (|2.1[) . Fig. [3] characterizes the in- 
fluence of the effective exchange anisotropy for h x ^ 
and h z = 0. When A e fr > 1 (A < 0) the only available 
ground state is 'W 2 ". When A cff < 1 (A > 0), "iV*" 
prevails at h x /D >> 1 and LL reigns at h x /D « 1. 
From Eq. (|4T6| it is clear that "A™" contains dimer- 
ization that is most noticeable near the phase boundary 
with LL, fading away as h x j D oo. As the easy-plane 
anisotropy gets stronger, the range of h x /D for which 



"N z " is the ground state becomes narrower. As a mat- 
ter of fact 'W 2 " disappears completely for A c h < A c , 
where the critical value A c depends on h x , D and h z . 
We shall revisit and reafirm this point in the next sub- 
section. For A c ff < A c , the phase boundary between LL 
and U N V " is independent of A and occurs at h x ~ 1.5D. 
This regime matches fairly well with the classical pre- 
dictions of Section III, where the critical field for the 
commensurate-incommensurate transition was found to 
be h x , c = ttDS = ttD/2. 

Fig.@]differs from Fig.|3]only quantitatively, but serves 
to highlight that 'W 2 " is more robust for larger values of 
D and h x even as h x /D is unchanged. This observation 
can be understood as follows. When A e ff < 1, the value 
of yc decreases during the first stage of RG. When D 
and h x are very small, 1$ is large and the prerequisite 
for flowing to 'W 2 " (yc(l<t>) > -y*(l<t>) at y a < 0) is less 
likely to be fulfilled. The larger D and h x are, the smaller 

and thus the shorter the decay of yc ; this improves the 
odds for a 'W 2 " ground state. 

Fig. \5\ shows that for a given anisotropy the phase 
boundaries between LL, "N y " and 'W 2 " are linear. It 
is straightforward to explain this behavior analytically. 
Since h z = we have 9 + = —tt and hence ?/a(0 = 0. 
Furthermore, we shall assume that yA is small in the first 
stage of the RG flow. This assumption is adequate only if 
sin 6~ ~ 0. Fortunately, this condition is satisfied when 
D » h x (9~ ~ 7r) or D « h x (8~ ~ 0), which arc 
the regions of interest when seeking u N yv or LL. Thus 
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TABLE I: Ground states for Eq. (|2.1[) . based on the value of the coupling constants at the end of the RG flow. Our calculations 
show that {£) = for any h x ,h z and D. However, this does not preclude a dimerized phase in Eq. (|2.1[) because the gapped 
phase labeled as "N v " contains a mixture of Neel-y correlations along with dimerization. Our RG analysis demonstrates that 
u N y " emerges when D « h x and A < 1. It follows from Eq. (|4.16|) that the dimerization component of this phase is enhanced 
as D gets larger. The nomenclature for gapped phases labeled as U N X " and "iV z " is motivated by the h z — case, for which 
"N x " ("JV Z ") means antiferromagnetic order with Neel vector along x (z). When h z ft the Neel vector lies in the xz plane 
for both phases. 
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FIG. 4: Phase diagram corresponding to Eq. f|2 . 1 j) : derived 
using non-Abelian bosonization. D = 0.01J and h z = 0; we 
vary h x and A. Comparing this plot with Fig. [3] it is clear 
that for a given D/h x , "A z " is more robust when D is larger. 
In other words, A c (~ 0.94 in this figure) decreases as D and 
h x increase. 



neglecting yA and ija from the onset, we are left with KT 
equations for yc and y a , the initial conditions being 



tfc(0) 



gbs 
2ttv 



1 



ffbs 

Attv 



1 



cos 9 



cos 9 



1 



(4.20) 



On one hand, the KT equations lead to "N y " provided 
that yc < and y a > (or alternatively if yc < y<r 
and y a < 0, although this is not satisfied at small A). 
After some algebra this amounts to requesting 8~ G 



(-2VA,2VA), for A << 1. Using D/h x = tan0 _ ~ 6~ , 
it follows that exists when < D/h x < \/A. 

This is why we have a linear phase boundary between 
"N v " and LL. Moreover, the slope of the correspond- 
ing line in Fig. [5] matches fairly well with VA. On the 
other hand, the KT equations predict LL if y a < and 
VC < Va < —yc- These inequalities may be reduced 
to cos8~ < 2 A — 1, where we have expanded for small 
A. This condition is satisfied for 9~ £ (n — 2^/\,tt) or 
equivalently D/h x > VA, which is another straight line. 



0.02 



0.015 



-x 0.01 



0.005 



0.02 




FIG. 5: Phase diagram with A ~ c(l — A e g) ~ 0.02 and 
h z = 0; derived using non-Abelian bosonization. "N y " pre- 
vails when D/h > y/\ and LL prevails when h/D < \T\. 
These results may be understood analytically, as discussed in 
the text. 



For completeness we include Figs. El [7] and HI which 
display typical RG flows for the coupling constants in 
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FIG. 6: (color online) Typical flow of the coupling constants 
in the LL phase. D = 0.01J, h x /D = 0.1, h z = and A ~ 
0.975. Va(1) = and j/b(£) — (not shown). A horizontal 
line indicates that a particular coupling constant has stopped 
flowing at because of rapid spatial oscillations, ~ 2, 1$ = 

CO. 



FIG. 8: (color online) Typical flow of the coupling constants 
in the "iV z " phase. D = 0.01J, h x /D = 1.5, h s = and 
A ~ 0.975. ju(Z) = 0(not shown). 1$ ~ 1.5, Z e = oo. Had 
we chosen a smaller A the strong coupling would have been 
reached at larger value of I. This is why the "JV Z " phase at 
A < 1 is particularly fragile to h z / 0. 
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FIG. 7: (color online) Typical flow of the coupling constants 
in the "N y " phase. D = 0.01J, h x /D = 5, h z = and 
A ~ 0.975. va{1) = 0(not shown). 1$ ~ 0.4, l e = oo. 



values of Z) and /i^.— This means that the correlation 
length of the antiferromagnetic states is £ > aoexp(10). 
When h z ^ 0, tg is finite and insofar as tg < exp(— 10)a ( ^ 1 
the flow of yc stops before it reaches the strong coupling 
limit. Specifically, the critical value for the field is 



(4.21) 



which wc plot in Fig. (fTTj) . All in all, this figure indicates 
that the field-induced long-range order ("iV z ") is hardly 
detectable experimentally when A < 1. For small values 
of h z , the typical flow diagrams for h z ^ are identical 
to those of Figs. [51 [7] and |S] except for the fact that now 
the plots must end at / = lg. As an example, consider the 
parameters of Fig. [5] with h z — 10~ 6 h x . It follows that 
lg ~ 16. The ground state is LL because |yc(16)| << 1, 
as shown in Fig. [5] 



each of the phases. 

Thus far we have taken the magnetic field to be 
completely perpendicular to the DM vector. Figs. [S] and 
1101 evidence that even a tiny h z can modify the phase 
diagram substantially, owing to the large correlation 
length of the ordered states near A = 1. Let us denote 
as l c the value of I at which \yc\ — 1. It follows from 
the analytical solution of the KT equations^ that 

lc = [tt - cos- 1 {-y a m)/yc{h))\l Wycih) 2 - y<r{l<f,) 2 
(for \yc{l<j>)\ > \vJU}\) or lc = 

co S h-\yS^)/\vc{h)\)/Vy^? ~ Vcih) 2 (for 
\yc(l<t>)\ < \y<r( l <t>)\)- We nnd l c > 10 for typical 



B. Abelian bosonization 

The main limitation of the non-Abelian bosoniza- 
tion method employed so far is that it is designed for 
SU(2) symmetric systems, thus becoming gradually un- 
reliable away from A e g = 1. In this section we at- 
tempt to overcome this limitation by turning to Abelian 
bosonization, which is better suited to handle exchange 
anisotropy. Nonetheless, the Abelian bosonization ap- 
proach has shortcomings of its own. Most notably, there 
is no obvious way to perform the chiral rotation that 
proved very helpful in the non-Abelian case. All in all, 
this subsection does not presents new results but rather 
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FIG. 9: Influence of h~ on the phase diagram of Eq. (|2,1[) ; 
derived using non-Abelian bosonization. h z = 10~ 6 h x and 
A ~ c(l — A c ff) ~ 0.02. The LL ground state is defined 
(somewhat arbitrarily) via |yc('e)| < 0.2. Even for a small 
2-component of the field the outcome is dramatically differ- 
ent from Fig. In this figure "N y " is no longer present (it 
resurfaces for smaller A e ff) and "7V Z " can be found only in a 
limited region of the parameter space (that of D ~ h x and 
relatively large D). The reason why U N V " is more fragile than 
l 'N z " is that the value of / at which it reaches strong coupling 
is larger. As we shall see in the next subsection, this trend 
reverses when the easy-plane anisotropy is stronger; in that 
case "N z " is more fragile than U N V " . 




FIG. 10: Influence of h z on the phase diagram of Eq. (|2,1[) ; 
derived using non-Abelian bosonization. h, = 10~°h x , D — 
0.01J. The LL ground state is defined (somewhat arbitrar- 
ily) via |f/c('e)l < 0.2. Comparing this figure with Fig.[4l it is 
clear that even a small h z brings about qualitative changes to 
the phase diagram. Remarkably, not only "iV z " is no longer 
the ground state for A < 1, but even for A > 1 there is a 
swath of parameter space for which LL prevails. For A c ff < 1 
the phase boundary between LL and "JV 9 " is pushed to larger 
values of h x . This is in qualitative agreement with Eq. (|3.1|l : 
the critical field for the commensurate-incommensurate tran- 
sition increases due to h z , such increase being less pronounced 
as A is made smaller. 




2e-06 



FIG. 11: Plot of Eq. (^2TT) for A = 1 and D ~ 0.1J. When 
h z > h Ct z "jV z " is no longer the ground state. /i ZjC is largest 
when h T ~ D. 



revisits from a different viewpoint the overall features of 
the phase diagram derived above. 

We begin by gauging away the DM interaction as ex- 
plained in Section I; this renormalizes the anisotropy pa- 
rameter A — > A c ff = A/-\/l + D 2 /J 2 and introduces 
spatial oscillations in the component of the magnetic field 
that is perpendicular to the DM vector. Next, by media- 
tion of the Jordan- Wigner transformation we express Sj 
in terms of spinless fermions ipj i Eq. (|2.1[) then turns into 
n=n +n hs + V, where 



Ho = J 



(4.22) 



is the non-interacting part, 



H hs = JA cS J2 UHo - I) (Vj+i^+i - I) (4-23) 



is the interacting part (A e g plays the role of interaction 
strength) , 



V 



- hz y: (frfi ^ 

h * (y]e laj e™ Efc <^ ^ k + h.c.) (4.24) 



collects the DM and Zccman terms and a = tan _1 (Z?/J). 
The low energy properties of this model may be captured 
by linearizing the dispersion of the spinless fermions 
around the Fermi points and going to the continuum 
limit: 

tpj/y/E ~ Re lkFX + Le- ikFX ; tp j+1 ~ ijj(x) + ad x ^{x), 

(4.25) 
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where kp = it /2a + h z /v is the Fermi momentum (half- 
filling at h z = 0). The right- and left- moving fermion 
fields arc bosonized via 



R = 



1 



1 



-i V47T0i 



(4.26) 



/2na v27ra 
where 4>R,l are chiral bosonic fields obeying 

i 

[4>Ri 0l] = J 

[4>R(L)(x),4>R(L){y)} = +(-)|sign(a: - y). (4.27) 

Symmetric and antisymmetric combinations of the chiral 
fields constitute the dual fields introduced in the previous 
subsection: 

§ = <t> R + 4>L ; 9 = -4>R + <t>L- (4-28) 
Using $ and the bosonized form of Eq. (|2.ip reads 



H = H Q + H hs + V 



Ho 

Hbs 

V 



2ttw ^ / . 2h z VK 

r Gcos V8ttA $ H — — x 



{2iray 



h 
na 



- cos ( V2ttK$ + — ^-^-x 




X cos 



In the derivation of V we have used 



(4.29) 




(-1)*/° + cos(V2tt#$) 

iW a , 

— cos(V27rif$), (4.30) 



neglecting rapidly oscillating terms and absorbing h z and 
D through a shift in $ and 0, respectively. In addi- 
tion, we have rescaled the bosonic fields as <& — > vA 7 ^? 
and — >• Q/y/~K, which enables us to write the non- 
interacting part Ho in the canonical form. K is the 
Luttinger parameter that differs from unity due to in- 
teractions (A c ff 7^0). = 2(1 — K), which played 
a central role in the preceding subsection, characterizes 
the interaction strength. The cos(v8tt^$) term orig- 
inates from Umklapp scattering events such as R'(x + 
a)L(x + a)L'(x)R(x). When h z ^ the fermionic sys- 
tem is away from half-filling, hence the onset of spatial 
oscillations. Likewise, D induces spatial oscillations in 
cos{\j '8ir I 'KQ) , which can be understood by carrying out 
the bosonization in a rotated frame.— 

The Umklapp perturbation has zero conformal spin 
and its scaling dimension is d = ( \/8nK) 2 /4n = 2K, 
which is irrelevant for K > 1 (A c g < 1). In contrast, V 



contains a product of two operators with nonzero confor- 
mal spin s = V2^K v / 2tt/K/2tt = 1. As a consequence, 
V generates new perturbations in the course of the RG 
flow, which have to be taken into account. These pertur- 
bations have zero conformal spin and may be derived as 
indicated in Ref. [13]; the outcome is 



n = - [(d x ef + a*) 2 ] 



4vz / , h z VK 

H — 7T- cos \/2-kK<$> + — — 
or \ v 



x cos 




K aVK' 



^-Gcos V&kK® + — ^ — j 

(2-Kay \ v 



2ttv 
' (2?ra) : 



-G cos 




8tt 2a 

— e + — ^x 



(4.31) 



where 



= hxa/Anv. The new perturbation gener- 
ated from V is cos( ^JR-k/KQ), with scaling dimensions 
d = 2/K. This perturbation is relevant at A g < 1, 
which combined with the fact that cos(V87rA$) is irrel- 
evant suggests that the ground state should be described 
by a pinned field. This guess is naive, partly because 
when 2{K - l/K) < 1 (namely K G (1, 1.28)) V must 
be considered together— with the spinless perturbations, 
which complicates the outcome. Moreover, there are the 
oscillatory phases that stop the flow of the coupling con- 
stants. As in the non-Abelian study, we elicit the ground 
state of H from a multiple-stage renormalization group 
analysis. 

In the first stage of RG the characteristic momenta are 
larger than h z /v and a/ a (hereafter we neglect factors of 
order one that multiply the oscillatory phases) and we 
can set h z = D = in Eq. i|43I]) . When h z = D = 
Eq. (|4.3ip describes two weakly coupled Luttinger liq- 
uids; the corresponding flow equations are knownil to 
be 



dz 
~dl 
dG 

~dl 

dG 
~Jl 
dK 

~dT 











-±): 


z 









— = 2(1 -K)G 



= T 

K 
~2 



K\z A 



1 



G 1 



K 



G 2 K 



(4.32) 



When z = and K > 1, G flows towards strong 
coupling and G flows to weak coupling. If G reaches 
~ 1(— 1), the ground state is described by an ordered 
(disordered <3>) such that cos(y/8njKQ) = 1(— 1). From 
Eq. (|4.30p . this implies that only the staggered compo- 
nent of S x (S v ) acquires a nonzero expectation value. 
Consequently, the ground state is "N x " ("A y "). Con- 
versely, when z — and K < 1, G flows towards strong 
coupling and G flows to weak coupling. If G reaches ~ 1 , 
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the ground state is characterized by an ordered $ (disor- 
dered 9) such that cos(V&ttK$) = 1. From Eq. (|430|) . 
this implies that only the staggered component of S z ac- 
quires a nonzero expectation value. Consequently, the 
ground state is "A z ". When z ^ the aforementioned 
trends are less clear, and a more careful analysis is re- 
quired. We note in passing that for a = z = Eq. (|4.31[) 
is the XYZ Thirring model, which as h z increases un- 
dergoes a series of phase transitions from a commensu- 
rate (gapped) spin density wave (SDW) to an incom- 
mensurate (gapless) SDW back to a commensurate SDW 
through a spin flopi 17 ' 20 

Eq. (|4.32[) is most reliable for A ff ~ and 
h x ,h z ,D « J, because under these conditions the 
initial values for z,G, G and y a are guaranteed to be 
small. In effect, z(0) oc h x /J, G(0) oc A e ff + (3z 2 and 
G(0) oc z 2 , where (3 is a constant that may be derived 
perturbativelyJ^ A(0) can be reliably determined by in- 
tegrating Eq. (|4.32[) backwards so that for h x = h z = 
one reproduces the well-established LL ground state with 
A'(oo)- 1 = 1 - cos-^Aeff)/^ and G(oo) = G(co) = 
z(oo) = 0. For weak magnetic fields, one may still use 
the same A(0) to a good approximation. When A ff ~ 0, 
K renormalizcs little and hence A(0) ~ 1 — 2A G ff/7r. 

In any event, we are most interested in accessing the 
strongly interacting regime A c ff < 1 so that we can 
make contact with the previous subsection. In partic- 
ular we wish to find out how robust the field-induced 
"iV z " ground state is when the easy-plane anisotropy in 
enhanced. Unfortunately, for A e ft < 1 the initial values 
for the coupling constants are uncertain. The underlying 
reason is that G(0) ~ 1, which renders Eq. (|4.32j) in- 
valid. A more sensible approach is to assume that there 
has been some prior renormalization group flow (with un- 
known flow equations), which starting from strong cou- 
pling has culminated in a relatively small value of G at 
some I = Iq. The rationale behind this assumption is 
that for A e ff < 1 the cos(v87rA$) term is irrelevant. 
Thereafter Eq. (|4.32|) determines the flow at I > lo, and 
we are left to guess the initial conditions of the cou- 
pling constants at I = Iq. We take z(Iq) oc h x /J and 
G(?o) oc z(Iq) 2 , with proportionality constants of order 
unity. On the other hand, we choose the value for K(Iq) 
by hand; this is tantamount to selecting an intermediate 
energy scale for I = lo, which corresponds to a lengthscale 
ao that is larger than the lattice constant a. Finally, we 
integrate Eq. ()4.32|) backwards to determine G(Iq) such 
that the ground state in absence of fields will reproduce 
G — > and K — > Aj n f. Since h Zl h x << J, presumably 
the value of G(lo) will be nearly independent of the mag- 
netic field. Overall, our choice of K(l ) is engineered in 
a way that Eq. (|4.32[) will reproduce the known ground 
states of a variety of limiting cases, without having to 
tune the initial values for the coupling constants. These 
limiting cases are 

(i) A eff < 1, h x = h z = 0, any D. 

In this case z{l) = G(l) = and we are left with flow 
equations for G and K. We find that the ground state is 



LL, in agreement with Bcthc ansatz calculations. 

(ii) A cff < 1, h z = D = and h x ^ 0. 

In this case G flows to strong coupling: G —> — co, 
G — > 0, K — > oo. This corresponds to the U N V " ground 
state, which is the expected answer as discussed in the 
previous subsection. 

(iii) A ef f <l,h x = 0, h z ^ 0, any D. 

Here z(l) — 0. In the first stage of RG (I < h = 
min(log(i>/ao/i z ), log(a/aoo;)) G is the relevant perturba- 
tion. However, because G(Iq) oc z(Iq) 2 and z(Iq) oc h x , we 
have G(l ) = and dG(l )/dl = 0. Therefore G(l) = 
and there is no possibility for a "N v " ground state. The 
effective RG equations for the first stage are thus 



— = — K 2 G 2 . (4.33) 
dl 2 K 1 

Since K(Iq) > 1, G decreases (it is irrelevant at A e ff < 1); 
so does K, but more slowly than G. Therefore G cannot 
reach the strong coupling limit either and moreover its 
flow stops at I = if. In sum, the ground state is LL, 
which agrees with Bcthe ansatz results. 

(iv) A cff = 1, D + 0, h x + and h z = 0. 

A e ff = 1 is the situation for which the non-Abelian 
bosonization discussed above is reliable. In this case z 
and G stop flowing at I = log(a/aoa), beyond which we 
can set z = G = in Eq. (|4.32[) and keep integrating 
the flow equations for K and G. We obtain G — > co and 
K —> —co regardless of the h x /D ratio. This corresponds 
to the "A z " ground state and is in agreement with the 
results derived in the previous subsection. 

Next we compute the more general phase diagram nu- 
merically. We begin by taking h z — 0. Fig. [12] displays 
the three phases that compete with each other when 
A ff < 1, and agrees roughly with Fig. U derived from 
non-Abelian bosonization. Below A c ff < A c , "A z " is no 
longer the ground state; instead, u N y " and LL are sta- 
bilized depending on the h x /D value. When D = and 
h x 7^ 0, the ground state is invariably U N V " . Likewise, 
when D =/= and h x = 0, the ground state is invari- 
ably LL. However, "A z " can be the ground state at 
D ^ ^ h x and A e ff < 1, even though cos(V87rA$) 
is irrelevant, because of the frustrating influence that D 
and h x have on each other. Similarly, a LL ground state 
may emerge at h x ^ due to the DM interaction, which 
may cut off the flow of G before it reaches the strong cou- 
pling regime. We use G(li) > 0.15 as the criterion that 
defines "A w " (note that G(l ) « 1). This is a some- 
what arbitrary choice that endows the phase boundaries 
with uncertainty; nevertheless it is motivated as an at- 
tempt to reproduce the phase diagrams of the previous 
subsection. 

Unlike in the previous subsection, in the present con- 
text the "A a " phase displays no trace of dimerization. 
In effect, here dimerization is associated with the order- 
ing of $ such that cos(V87rA'$) = -1; yet in the il N v " 
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FIG. 12: Influence of exchange anisotropy on the phase 
diagram of Eq. (|2.1I) ; derived using Abelian bosonization. 
D — 0.01J, h z — 0. The initial energy scale for the RG 
flow was chosen via K(lo) = l.lA(oo). For A e ff < 0.97 only 
u N y " and LL phases can be found. This plot agrees roughly 
with Fig. 21 which was derived using non- Abelian bosoniza- 
tion. Disagreements between the figures are most noticeable 
on the shape of phase boundaries. 




0.005 0.01 0.015 0.02 0.025 
D/J 



FIG. 13: Phase diagram for h z — 0, A = 0.93; derived using 
Abelian bosonization. The initial energy scale for the RG flow 
was chosen via K(lo) = l.lA'(oo). For such A, "A"" is absent 
for D < 0.015 J (see also Fig. I12p : it reappears at larger D 
values. This demonstrates that A c gets smaller as D and h x 
increase. 



phase <E> is disordered because is pinned. The fact that 
non- Abelian bosonization is able capture the coexistence 
of antiferromagnetism and dimerization is more a merit 
of the chiral rotation than an intrinsic flaw of the Abelian 
bosonization. 

Fig. [13] sheds light on the parameter space for which 
l W z " constitutes the ground state. For given A e g < 1 
and h x /D, N z is more robust at larger D (or h x ). In 
other words, A c decreases as D and h x increase and their 




FIG. 14: Phase diagram in the neighborhood of the non- 
interacting limit (A = 0.1), where the Abelian bosoniza- 
tion scheme utilized is most reliable. "N y " is defined via 
G(h) > 0.15; other choices would lead to a quantitative 
change in the slope of the phase boundary between u N y " and 
LL. Regardless, there is no remnants of the "A z " phase. 



ratio is kept of order one; this is in agreement with the 
results derived in the previous subsection. The reason 
behind this trend is that G and G compete which each 
other, the latter being the most relevant perturbation at 
A c ff < 1. The larger D is the less chance G has to reach 
strong coupling, thereby increasing the likelihood for G 
to prevail. 

When h z ^ 0, the flow of G is interrupted. In qual- 
itative agreement with Fig. [5J we find (not shown) that 
the parameter range for "N*" narrows significantly even 
for a small value of h z . However, we remark that in the 
present formalism h z hampers "N z " but not "N v " be- 
cause h z does not induce oscillations in the factor that 
multiplies G. This is in qualitative disagreement with 
the outcome of non- Abelian bosonization, where h z ^ 
reduces the likelihood of both "A™" and U N Z " (because 
both phases are linked to the same coupling constant 
yc, cf. Sec. IVa). The classical study of Section III 
rules in favor of the non- Abelian result by anticipating 
an increase of the critical field for the spin-flop transition 
when h z ^ 0. 

As mentioned above, the main results of this subsec- 
tion suffer from uncertainties in the initial conditions for 
the coupling constant. This problem is remedied in the 
weakly interacting limit (A ft ~ 0) for which the initial 
values of the coupling constants are small and known. 
Fig- HH demonstrates that for large easy-plane anisotropy 
(A = 0.1) only a N yn and LL can be the ground states. 
This, in conjunction with the non- Abelian study of the 
SU(2)-symmetric point, ratifies that there exists a criti- 
cal value of A below which "N z " disappears. While our 
calculations indicate that A c is close to one, numerical 
density-matrix RG studies might be desirable to ascer- 
tain its precise value, as well as to corroborate the co- 
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existence of dimerization and antiferromagnetism in the 
ll N v " phase. 



V. APPLICATION TO SUPERCONDUCTIVITY 
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Thus far we have discussed the ground states of 
Eq. (|2.1[) in the context of one-dimensional quantum spin 
chains. However, there appear to be few experimental 
studies on one-dimensional antiferromagnets with uni- 
form DM interactions. In this section we demonstrate 
that Eq. (|2.1j) also models superconducting nanostruc- 
tures that might be realizable in experiments. 

First, consider a onc-dimcnsional array of Joscphson 
junctions^ separated from one another by a distance a. 
Its Hamiltonian is 



H 



1. 1 



2 ^ 



n i G i,j n j - E J^2 cos (Xj ~ Xj+l), (5-1) 



where rii denotes the number of Cooper pairs in the i-th 
superconducting island, \i is the U(l) superconducting 
angle for the i-th island (canonically conjugate to n^), 
Cjj is the capacitance matrix that models the repul- 
sive Coulomb interactions between the Cooper pairs, and 
Ej is the Josephson coupling energy. For conventional 
Joscphson junctions Ej > 0, while for 7r-j unctions^ 
Ej < 0. We neglect dissipativc processes (e.g. quasi- 
particle tunneling), which are relatively unimportant at 
low temperatures. We are interested in small supercon- 
ducting grains^ where the onsite Coulomb interaction is 
strong, i.e. e 2 C~l >> Ej,T. Accordingly the super- 
conducting islands are in the Coulomb blockade regime 
and the large electrostatic energy cost for changing the 
number of Cooper pairs on each island drives the JJ ar- 
ray to an insulating regime. We consider the particu- 
lar case in which rii can acquire only two possible val- 
ues; the practical implementation of this scenario may 
require tuning the chemical potential of the Cooper pairs 
via a gate voltage. The two possible values of charge 
define a pscudospin degree of freedom for each island, 
which enables the mapping of Eq. (|5.1[) into a pscudospin 
Hamiltoniani2£>2& 



Ujj = -EjJ2(S*S* +1 + Sp» +1 ) 

3 



E, 



E 

3 



QZ QZ 



(5.2) 



Sj = rij is the number operator for Cooper pairs, 
and E c Sj +1 Sj +1 describes the intcrgrain Coulomb re- 
pulsion (E c > 0). We have assumed screened Coulomb 
interactions rl whereby C~J » C~~J +l and C~l +n = 
for n > 2. This requires that the self-capacitance of 
the superconducting island Co be larger than the junc- 
tion capacitance C. Co and C arc defined via Ci.j ~ 
(C + 2C)S itj - C{S ilj+1 + Sij-i). S* and are the' real 
and imaginary part of the superconducting pair operator 



FIG. 15: Superconducting analogue of Eq. (|2.1[) : a one- 
dimensional array of small superconducting islands (in black) 
separated by Josephson junctions (crosses), located in close 
proximity to a bulk superconductor (shaded area). The ap- 
plied magnetic field (perpendicular to the page) leads to an 
effective DM interaction in the JJ array. The Josephson cou- 
pling between the bulk superconductor and the array plays 
the role of a XY magnetic field in pseudospin space; gate 
voltages are Z magnetic fields in pseudospin space. The pseu- 
dospin anisotropy is defined by the disparity between the 
Josephson coupling and the capacitive energy of the junc- 
tions. 



for the j-th grain, and Ej(SjSj +1 + SjSj +1 ) is associated 
with the tunneling of Cooper pairs between neighboring 
grains. h z is a pseudospin magnetic field that describes 
the deviation of the chemical potential from the middle 
point between the electrostatic energies of the two charge 
states. Unlike in quantum antiferromagnets, in JJ arrays 
Ej/E c may be tuned in situ^ 

Second, let us place a large superconductor parallel to 
the array of junctions (see Fig. [15]), separated by a dis- 
tance d. If the material placed between the JJ array and 
the bulk superconductor is a normal metal, a Josephson 
coupling will ensue as long as d < £j\r, where £jv is the 
coherence length of the normal metal. In magnetic lan- 
guage, the influence of the large superconductor is equiv- 
alent to that of an external magnetic field oriented in the 
xy plane: 



H = U 



j j 



E^ 



(5.3) 



where our gauge choice is determined by h y = 0. h x 
is proportional to the mean field order parameter of the 
large superconductor. Unlike in the JJ array, we neglect 
phase fluctuations in the bulk superconductor. Eq. (|5.3[) 
is a XXZ model with a uniform pseudo-magnetic field. 

Third, we add the ingredient which will result in a 
Dzyaloshinskii-Moriya interaction in pseudospin spaced 
Let us apply a uniform magnetic field B = Bz (see 
Fig. IT5|) . The vector potential A associated with the 
magnetic field twists the superconducting angle, so that 

E (*"stsr + i + h - c -) + & E s i s Ui 



n 



ii-Y,(<" !,s ; • '-••) ^Es 



(5.4) 



where 8 = (2tt/$ ) f- +1 Mv = 0) ■ dx and 



(2tt/3> ) J j+1 Mv = d ) ■ dx - $o = h/2e = 2x 10~ 15 Wb 
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is the flux quantum and we have taken y = at the loca- 
tion of the J J array. Eq. (|5.4[) assumes that in spite of the 
vector potential the superconducting phase is spatially 
homogeneous within each island. This is a reasonable ap- 
proximation insofar as the magnetic flux threading the is- 
land is small. For convenience we use A = —Byx, which 
results in 



H 



E 



(a) 



t f / / 



(b) 



i i 



® 



111! 

( C ) ! fYo E fY 

W W 



(5.5) 



with 



<f> R 

a = -27T ; $ B = Bda. 

$ 



(5.6) 



(f>s is the magnetic flux penetrating a rectangle formed 
by a Josephson junction, perpendicular lines to the bulk 
superconductor and the edge of the bulk superconductor. 
Our choice of the vector potential corresponded to a 
spatially uniform magnetic field; nevertheless, Eq. ([5. 5ft 
is valid more generally. If Ej < (7r-junctions), 
Eq. (|5.5j) completes the mapping into Eq. (|2.2p . In 
contrast, if Ej > (conventional junctions) we need to 
make an additional pscudospin rotation by an angle tt 
for every other site: ST — > exp(iirj)Sj~ . The resulting 
correspondence between the antiferromagnetic and the 
superconducting models can be summarized as follows: 

Antifcrromagnct Conventional JJ Array ir - 33 Array 



V J 2 + D 2 
JA 

a = tan~ 1 {D/J) 



Ej 
E c 
a + 7r 



-Ej 
E c 
a 

hrr. hv 



The phase diagrams calculated in the previous sections 
are valid for D,h x ,h z << J. Therefore, those results 
may be transferred directly to the case of conventional 
junctions only when 4>b ~ $0/2 and h Xi h z << Ej. 
On the other hand, for 7r-junctions our RG analysis has 
access to $b ~ and h x ,h z « Ej. 

Fig. [16] illustrates the physical meaning of the mag- 
netic phase diagram in the present context. "N zv corre- 
sponds to an insulating charge density wave phase, where 
the number of Cooper pairs oscillates from one island to 
another. "N y " is a vortex phased where there are circu- 
lating currents with alternating chirality flowing between 
the JJ array and the bulk superconductor. In addition, 
"N v " contains some dimerization: 



(5.7) 



(S+Sr +1 ) = a + (-iy/3 

(S*S? +1 ) = j+(-iyr h 

where a,(3,j,r) are constants. The first line of Eq. ([57 
implies that the effective Josephson coupling between 
sites 2j and 2j — 1 is larger than that between sites 2j and 



(d) 



FIG. 16: Analogues of the magnetic phases in conventional 
Josephson junction arrays - classical representation for h z = 
0. Dashed arrows represent the direction of the order param- 
eter (aj) in the bulk superconductor. The gauge choice is 
determined by Eq. H5.5|) . Solid arrows portray the order pa- 
rameter of the superconducting islands in the one-dimensional 
JJ array. The vector potential twists the direction of the pseu- 
dospins in the XY plane, (a) Ferromagnetic phase. This is 
the classical ground state for small magnetic fields ($_b — ► 0), 
a regime in which the field theoretical results of the present 
paper do not apply. The spatial gradient of the order pa- 
rameter implies Meissner currents flowing at the edge of the 
bulk superconductor, (b) "iV z " phase at $s — $o/2 (for 
simplicity we have plotted the order parameter of the bulk 
superconductor as though a = n). This charge density wave 
state is dominant at E c > Ej and may also arise at E c < Ej 
provided that |$s/$o — 1/2 j ~ h x /Ej (although the experi- 
mental detection in the latter regime is highly unlikely). The 
pseudospins of the array point along Z with alternating sign, 
(c) "N v " phase at $s ~ $0/2 (for simplicity we have plotted 
the order parameter of the bulk superconductor as though 
a = 7r). There are circulating supercurrents of alternating 
chirality and oscillating magnitude across the interface be- 
tween the bulk superconductor and the ID array, (d) LL 
phase at $s ~ $o/2 (for simplicity we have plotted the or- 
der parameter of the bulk superconductor as though a = ir). 
When $_b — $o/2 = the array ignores the twist in the order 
parameter of the bulk superconductor and adopts a ferromag- 
netic configuration (not shown) along an arbitrary direction 
in the XY plane, much as though the bulk superconductor 
did not exist. When ^= |$s — <&o/2| < h x /Ej the classical 
configuration is a soliton lattice (only one soliton is shown in 
the figure). 



2j + 1. In other words, the magnitude of the circulating 
currents oscillates from one "plaquette" to another and 
is larger for one chirality than for the opposite chirality. 
The second line of Eq. (|5.7[) means that the magnitude 
of the effective junction capacitance oscillates from one 
plaquette to another. In other words, if island 2j has 
zero Cooper pairs then site 2j — 1 is more likely to have 
one Cooper pair than site 2j + 1. Finally, LL is the 
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gaplcss ground state in which the superconducting an- 
gles of the junctions form a soliton lattice. Throughout 
the foregoing discussion we have generally ignored the 
back action of the ID array on the bulk superconductor, 
and in particular we have neglected the magnetic fields 
generated by the alternating currents in the a N v " phase. 
This approximation is safest when the effective London 
penetration depth is larger than the size of the array^i 

Let us discuss the phase diagram of a conventional ar- 
ray with E c < Ej, which in the magnetic problem corre- 
sponds to A cf j < 1. When $ B ~ 0, the order parameter 
of the JJ array is aligned ferromagnetically with that 
of the bulk superconductor. Although this result does 
not follow from the field theoretical calculations of the 
present paper, its counterpart in quantum spin chains 
is well-established. In effect, for a = and h z = 
the superconducting model becomes equivalent to a spin 
1/2 antiferromagnet with staggered DM interaction in 
a transverse magnetic field, which was first studied in 
Ref. @. The ground state in this case was found to be 
antiferromagnetic with the Neel vector aligned with the 
external field. In pseudospin language this translates into 
the aforementioned ferromagnetic ground state. 

As the magnetic field increases the pscudospins of the 
JJ array are increasingly twisted, with the concomitant 
loss of exchange energy. For $b ~ ^ci/2, it is no longer 
optimal to have a ferromagnetic alignment between the 
array and the bulk superconducting angles. Instead, 
the ground state is "iV"" (if |$ B /$ - 1/2| < h x /Ej) 
or LL (if l^s/f&o - 1/2] > h x /Ej). Furthermore, if 
E c /Ej G (A c ,l), ll N z " emerges as the ground state at 
|$s/$o - 1/2| ~ h x /Ej. Finding "iV z " at E c < Ej is 
counterintuitive because it means that the proximity cou- 
pling from the bulk superconductor drives the array into 
an insulating ground state. By tuning h x or the applied 
(real) magnetic field, one may induce transitions between 
the three phases. In particular, if E c /Ej < A c there is 
a commensurate-incommensurate transition between LL 
and 11 N y " . The critical value of the proximity coupling 
for this transition is 

/i XiC ~ l(a + ir)Ej. (5.8) 

This critical field changes in presence of a gate voltage 
(h z ) as indicated by Eq. (|3.ip . 

The experimental detection of the aforementioned 
phases requires SQUID measurements^ which would 
target the circulating currents of "N v " ground states, 
as well as measurements of the critical current of the 
array^ which would be exponentially suppressed with 
the length of the array in the "7V Z " ground stated but 
not in the LL phase. We briefly comment on a number 
of additional experimental requirements: 

(i) The Josephson coupling between the islands in the 
array must be stronger than the coupling between the 
islands and the bulk superconductor, because our field 
theoretical results apply for h x << Ej. Moreover, our 
model applies for short-ranged Coulomb interactions in 
the array, i.e. for C < Co- 



(ii) The temperature of the system must be smaller 
than the gaps in the U N V " and "iV z " ground states: 
T << h x « Ej,E c . For typical values of the Joseph- 
son coupling (Ej < IK for low-temperature supercon- 
ductors) this requirement is most pressing at E c ~ Ej, 
where the energy gap associated with "iV z " is only ~ 
0.077£j exp(- 10) ~ 5(iK (recall Section IVa). As E c 
and Ej are made dissimilar the gaps may increase to 
0(h x ) < lOmK . 

(iii) The array must be long enough so as to reduce the 
quantum tunneling between degenerate ground states. In 
effect, "N y " and U N Z " each break a Z 2 symmetry: there 
are two degenerate "iV 2 " or "N y " phases that differ from 
each other only by a translation of a lattice constant. 

(iv) The superconducting grains must be small enough 
to justify the pseudospin 1/2 approximation. At the 
same time the area of a plaquette (a x d) should be 
large enough to enclose a flux $b = $o/2 using mag- 
netic fields that arc smaller than the critical field of 
the superconducting islands (note that this concern does 
not apply to arrays of 7r-junctions, for which we require 
$b ~ 0). The "inverted-T" shape of the islands de- 
picted in Fig. ITSI could help satisfy both conditions. For 
aluminum, a x d > 0.1 fim 2 would ensure that the applied 
field remains below the critical field^ 



VI. SUMMARY AND CONCLUSIONS 

We have evaluated the zero-temperature phase dia- 
gram of an antiferromagnetic spin 1/2 chain in presence 
of uniform Dzyaloshinskii-Moriya interactions, symmet- 
ric exchange anisotropy and arbitrarily oriented magnetic 
fields. We have used non-Abelian as well as Abelian 
bosonization, and have generally found qualitative agree- 
ment between the two schemes. When the two diverge, 
the former approach proves to be more reliable. Our 
calculations predict the emergence of three competing 
phases for spin chains with easy-plane anisotropy. One of 
them (phase (i)) is an antiferromagnet with its Neel vec- 
tor along the direction of the DM vector. This phase was 
introduced in previous workp£ whose scope was limited 
by the assumption of isotropic symmetric exchange. Our 
results indicate that phase (i) is unstable under weak-to- 
moderate easy-plane anisotropy: we have estimated the 
critical value of the anisotropy beyond which it disap- 
pears. This value is sensitive to the magnitude of the DM 
interaction, as well as to the magnitude and direction of 
the applied magnetic field. The two new ground states 
that occur as a consequence of symmmetric easy-plane 
exchange anisotropy are (ii) a dimcrized antiferromagnet 
with Neel vector perpendicular to both the DM vector 
and the magnetic field, (iii) a gapless Luttinger liquid, 
whose classical counterpart is a soliton lattice. Phase (ii) 
arises when the DM interaction is weak compared to the 
magnetic field component transverse to the DM vector; 
phase (iii) ensues in the opposite regime. Phase (i) may 
then be understood as an outcome of the frustration be- 
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tween competing phases (ii) and (iii); indeed it is most 
likely to emerge when the DM interaction is neither large 
nor small compared to the transverse magnetic field com- 
ponent. It would be interesting to verify and refine these 
predictions using the numerical density-matrix renormal- 
ization group method. 

Motivated in part by the scarcity of experiments 
on one-dimensional antifcrromagncts with uniform 
Dzyaloshinskii-Moriya interaction, we have searched for 
alternative systems where our calculations may be ex- 
perimentally tested. Thus we have mapped the origi- 
nal magnetic problem into a mathematically equivalent 
superconducting problem involving a one-dimensional ar- 
ray of Josephson junctions (either conventional or 7r-type 
junctions) in close proximity to a bulk superconductor. 
An applied perpendicular magnetic field plays the role 



of a uniform DM interaction. We have discussed the 
physical meaning of the magnetic phases in the super- 
conducting context, including that of dimcrization. The 
exquisite tunability of Josephson junction parameters in 
one-dimensional arrays may provide an interesting av- 
enue to probe and replicate the influence of DM inter- 
actions and magnetic fields in one-dimensional quantum 
antifcrromagnets with symmetric exchange anisotropy. 
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Appendix A: Influence of h z on the Classical Commensurate-Incommensurate Transition 



The objective of this Appendix is to determine how the classical soliton lattice (Sec. Ill) is modified when a 
magnetic field is applied along the direction of the DM vector. Let us parametrize the classical spin at site j as 
Sj = S(s'm9j cos fa , sin Oj sin fa , cos 9 3 ) . Then Eq. (|2.1[) can be rewritten as 

% = — JS 2 [sinflj sin#j + i cos(0j+i — fa — a) — A c ff cosdj cos#j+i] — h x S^^(— 1) J smdj cos fa — h z S cos 6> 3 ■, 



where 



7TJ, J = 3\f\ + D 2 / J 2 and A e ff = A/^/l + D 2 / J 2 . Taking advantage of the fact that fa+i 

1 , , 



(Al) 



is small for each j , we write 



sin 9j sin 9j+i 1 



a) I — A c ff cos 9j cos 9j+i 



-h x S'^^(—iy sin9j cos fa — /i 2 5^cos( 



(A2) 



Let us define 



fa = a(j) + (-iyb(j) 

o, = c(j) + (-iyd(j), 



(A3) 

where a, b, c, d arc functions that vary slowly along the soliton. 

Keeping only the non-alternating terms and making the continuum approximation we arrive at the following ex- 
pression for the Hamiltonian density h ("H = / dxh): 



-JS 2 ^[l-P 2 -d 
b 2 



-2b 2 

2 I dx 



- A eff (/3 2 - d 2 ) 



pd\l 



-h x S 



~h z S(3 1 



b sin a 



d 2 



where j3 = cose. In addition we have used smb ~ b, cos b ~ 1 — b 2 /2, sind ~ d, sine ~ 1 — ft 2 /2 and so on. 
Now we determine the optimal value for the functions a, b, c, d. The value of b that minimizes Hi is 

d 2 f3 2 " 



h x S sin a 
4JS(1 - P 2 - d 2 ) + h x $d cos a 
h r sin a 



1 



4J5 



B 2 d 2 h x Bd cos a 

1 + — — I — = 

2 2 4JS 



(A4) 



(A5) 
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Substituting this expression back in the Hamiltonian we obtain 



- JS 2 {(1 - /3 2 - d 2 ) 



1 1 ( da 

~ 2 \dx ~ a 



A cS {]3 2 - d z ) ) - h x S/3dcosa~ h z S(3 1 



h x sm a 
8J 



Next we optimize d. For A e ff < 1, Ih is minimized for d given by 



hx^cosa 
2JS(1 - A cff ) 



2(1 -A eff ) Vda; 



da 



— a \ — 



2JS{1 - Aeff) 



(A6) 



(A7) 



which indicates that the staggered component of 9 is nonzero only near the core of the soliton (far from the core cos a : 
as shown in Ref. (HI). For simplicity in the above approximation we assumed that |A c ff — 1| >> h 2 /J 2 S 2 , (da/dx- 
a) 2 . Substituting the expression for d in the Hamiltonian we get 



lh~ -JS 2 { (1-/3 2 ) 



1 



1 / da 



2 \dx 

Next we optimize /3. For A c ff < 1, Hi is minimized for /3 given by 



Aeff/3 2 } - h z S(3 



2J5(1 + A cff ) 



1 



1 



da 



h 2 sin 2 a h 2 (3 2 cos 2 a 
8J^ 4J(1 - A) ' 



h 2 cos 2 a 



2(1 + Aeff) \dx J AJ 2 S 2 (1- A 2 ff ) 



(A8) 



(A9) 



Plugging this expression back in Eq. IA8I and following with some algebra we get the sine-Gordon Hamiltonian with 
effective parameters: 



Hi = Hi 



'off-: 



o 



2 V dx 



2 h 2 cS sin 2 i 



-CM - 



8J e ff 



(A10) 



where 



[h = -JS 2 - 

Jcff = J 1 



h 2 



KK 



4J(1 + A cff ) 16J35 2 (l + A c ff)(l-A 2 ff ) 
hi 



4J 2 ,S 2 (1 + A e ff) 2 



h c « = h 3 



1 



h 2 



2J 2 5 2 (1 + A c ff)(l-A 2 ff ) 



1 - 



hi 



4J 2 5 2 (1 + A cff ) 2 



h 2 



3- Aeff 



AJ 2 S 2 (l + A off ) 2 (l- Aeff) 



(All) 



Following the same procedure as in Ref. [1J], the critical field for the commensurate- incommensurate transition is 
given by 



h e g = TraJeff, 



which yields 



hx,c 



3- Aeff 



8J 2 S 2 (1 + A cff ) 2 (l- Aeff) 

After some quick algebra this results in 



naJS 



u 2 s 2 (i + A eS y 



h x ,c = TTCtJS 



1 



8J 2 5 2 (1-A 2 ff ) 



(A12) 



(A13) 



(A14) 



which is precisely Eq. (13. 1|) . 



19 



Appendix B: Case Studies: Simple Regions of the Quantum Phase Diagram 



The objective of this Appendix is to verify the consistency of Eqs. (|4.11|) , (|4.13[) and Table I for a variety of cases 
in which Bethe ansatz solutions are available. As a byproduct we derive an expression for the constant c defined 
through A xc = c(l — A), and comment on our choice for the RG cutoff energy scale. 

(i) XXZ model with h x = h z = and D ^ 0. 

In this case 8 R = 0, 9 L = —n and thus ua{0) = 2/a(0) = 0. It follows that Da{1) = VaQ) = 0- Accordingly the 
pertinent RG equations become 

dyx dy y dy z 

-[r = y*yv->^r = y*y*->^r = y*yv> ( B1 ) 

with initial conditions given by y x (0) = —y y (0) — g\, s /(2-Kv) and y z (0) — </bs(l + X)/(2irv). y x (0) = ~y y (0) implies 
VxQ) = ~Uy(l) an d thus there is no instability towards "e" or "F". Moreover because y x (l) = ~VyiJ) we are left 
with Eq. (|4.14p for yc and y a , which applies when I < 1$ as well as when I > 1$. The analytical solutions of the 
Kostcrlitz-Thouless equations dictate that when A > (easy-plane anisotropy) the system flows to the gapless LL 
phase (yc — > 0), whereas when A < (easy-axis anisotropy) the system evolves to il N z " (yc — > oo). We reiterate 
that the 0(D 2 / J 2 ) term in the bosonized form of V (cf. Eq. (|4.3[) ) is crucial in order to get the correct answer for 
the case in which A = 1 and D ^ 0. In particular, for A > Bethe ansatz calculations prove that yc — > and 
y a — > 2(1 — A'inf), where = 1 — cos _1 (A e ff)/7r. With this in mind we evaluate the value of the constant c which 
enters the definition of A. We start by recognizing that Eq. (|B1[) implies y a (Q) 2 — yc(0) 2 = 2/er(oo) 2 — yc(oo) 2 with 



2/c(0) = 5bs/27rw and y a (0) = -y c (0)(l + A). This results in y a (<x) ~ -V2\g hs /2-Kv ~ -\/2c\/l - A cff .g bs /27ru, 
where we have neglected 0(X 2 ) and 0((1 — A)D 2 /J 2 ) terms. Comparing this with the Bethe ansatz prediction we 
obtain 

2 27nA 2 



c= ~ 7.66. (B2) 

V 71 " 9bs J 

(ii) XXZ model with h x = D = and h z + 0. 

In this case Or = 9^ = 0. Namely, the "chiral" rotation is simply the identity. In this case too 2m(/) = ju(0 = 
and the RG equations are given by Eq. (|B1[1 . The initial conditions are y x (0) = y y (0) = — 5bs/(27rw) and y 2 (0) = 
y y (0)(l + A). It follows that yc(l) = and Eq. (|Bip turns into KT equations for y B and y„. For A < the system 
flows to — > — oo, which in the original coordinates corresponds to ,l N z " (recall Eq. I4.16|) . ys reaches strong 

coupling (yB = —1) when I = l c = (2ttv / gh s )(^ — cos _1 (l + A)) / y— 2A for A < 0. However, h z interrupts the flow of 
yB at I = Icj, (note that tg = 0) and thus "iV z " is the ground state only when 1$ > l c ; for 1$ < l c the ground state is 
LL. Reaching strong coupling requires h z < h c where 



h c oc To exp 



2~KV 



9bs v / 2c(A - 1) 



(B3) 



is the critical field defined through 1$ = l c . The prefactor in Eq. (|B3[) is somewhat arbitrary because it depends on the 
precise value of J/b(Z) for which one decides that "strong coupling" has been reached. The critical field is also known 
from Bethe ansatz calculations which dictate h c / J oc cxp[— 7r 2 /2 y/2(A — 1)] for A > 1. Matching the exponent 
of this expression with that of Eq. (|B3j) yields c = 7.66, in agreement with Eq. (|B2j) . On the other hand, for A > 
Eq. (|B1[) flows to a weak-coupling regime regardless of h z (provided that h z << J). Once again this LL phase agrees 
with Bethe ansatz predictions. 

(iii) XXZ model with h z = D = and h x ^ 0. 

In this case 6r = 9l = —ir/2 and thus = Va(J) = 0. The flow equations are once again given by Eq. (|B1[) , with 

y y (0) = y z (0) — —gbs/(2Trv) and y x (0) = — <?bs(l + ^)/(2ttw). These RG equations must be replaced by Eq. (|4.14[) at 
I > Itp. We find that for A > the ground state is "A™" and while for A < the system flows to "iV z ". These results 
are in concordance with the classical considerations of Section I and agree with independent quantum mechanical 
calculations i^i We note in passing that there is no LL phase in the XXZ model with a uniform transverse field (i.e. 
D = in Eq. (j2~2i y 

(iv) XXX model with D = h x = and h z ^ 0. 

In this case y x (0) = y y (0) = y z (0) = — <?bs/27rw and the stage I flow equations reduce to dyi/dl = yf (for i = x,y, z), 
whose solution is yi(l) = —(gb s /2w)/(l + lghs/2-Kv). All couplings stop renormalizing at 1^ = log(v/aah z ). The 
Luttinger parameter K = 1 — \ya(l<j>) is then given by 
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Using To = v/ao = 0.077 J and g\, s = 0.23 x (2irv) we get ho = 5. 95 J. This is slightly smaller than ho = J\ZSn 3 /e = 
9.55 J, obtained from solving the Bethe ansatz equations.— However, this discrepancy is masked by the fact that the 
precise value of 1$ is uncertain (we could have used 1$ = log(Av j aoh z ) , where A is any constant of order one). 
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